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ABSTRACT 


We study nonlinear waves in a prominence foot using 2.5D MHD model moti¬ 
vated by recent high-resolution observations with Hinode/SOT in Ca II emission 
of a prominence on October 10, 2012 showing highly dynamic small-scale motions 
in the prominence material. Observations of Ha intensities and of Doppler shifts 
show similar propagating fluctuations. However the optically thick nature of 
the emission lines inhibits unique quantitative interpretation in terms of density. 
Nevertheless, we find evidence of nonlinear wave activity in the prominence foot 
by examining the relative magnitude of the fluctuation intensity (51 // ~ 5n/n). 
The waves are evident as significant density fluctuations that vary with height, 
and apparently travel upward from the chromosphere into the prominence ma¬ 
terial with quasi-periodic fluctuations with typical period in the range of 5-11 
minutes, and wavelengths ~<2000 km. Recent Doppler shift observations show 
the transverse displacement of the propagating waves. The magnetic field was 
measured with THEMIS instrument and was found to be 5-14 G. For the typ¬ 
ical prominence density the corresponding fast magnetosonic speed is ~20 km 
s^ 1 , in qualitative agreement with the propagation speed of the detected waves. 
The 2.5D MHD numerical model is constrained with the typical parameters of 
the prominence waves seen in observations. Our numerical results reproduce the 
nonlinear fast magnetosonic waves and provide strong support for the presence 
of these waves in the prominence foot. We also explore gravitational MHD os¬ 
cillations of the heavy prominence foot material supported by dipped magnetic 
field structure. 
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Introduction 


Large scale oscillations in prominences and filaments triggered by flares have been 


observed in Ha in the past using ground-based telescopes (e.g., Ramsey & Smith 1965 


Hyder 

1966; 

Thompson & Schmieder 

1991 

) as well as in He I ( 

Yi & Engvold 

1991) 


Observations of oscillations in prominences are important for prominence seismology 


e.g., Oliver & Ballester 2002 Oliver 2009 Arregui et al. 2012; Heinzel et al. 2014 


Ballester 2015), and due to the possible role of MHD waves in heating of the prominence 


material (Ofman & Mouradian 1996| |Ofman et al.| 1998). Recently with the advent 
of high resolution observations (both ground-based and space-based), detailed resolved 
observations of small-scale oscillations and waves in prominence threads became possible 
the review by Lin 2011]). In particular, Hinode/Solar Optical Telescope (SOT) (Kosugi 


see 


et al. 2007; Tsuneta et al. 2008) has provided high-resolution and cadence observations for 


transverse waves in prominences (Okamoto et al. 2007). However, the observations are often 


interpreted in terms of linear MHD waves using linearized models. The linear studies of 


prominence oscillations in terms of slab models were carried out in the past (e.g. Joarder & 


Roberts 

1992ajb, 

1993 

) as well as recently (e.g. 

Schmieder et al. 

2013 

Heinzel et al. 

2014) 


The magneto-acoustic gravity models in prominences have also been in the linear regime 


e.g., Oliver et al. 1992, 1993). While linear models may provide an adequate description 


of small amplitude waves in many observations, in some instances nonlinear effects can not 


be ignored, in particularly in large amplitude oscillations (Tripathi et al. 2009). Only few 


studies are devoted to theoretical models of nonlinear MHD waves in prominences (e.g., 
Chin et al. 2010| Terradas et al. 2013 Takahashi et al. 2015). 


Recently, Schmieder et al. (2013) reported detailed observations of propagating 


wave-like features in a quiescent prominence pillar (foot) observed by Solar Dynamics 


Obesrvatory (SDO) Atmospehric Imaging Assembley (AIA) (Pesncll et al. 2012 Lemen 
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et al. 2012) in EUV on 2012 October 10 with the details in high-resolution using the Hinode 
SOT Ca II and Ha lines, as well as ground based observations at Sac Peak and THEMIS 
(Telescope Hcliographique pour V Etude du Magnetisme et des Instability Solaires). The 
THEMIS observations utilized the MTR (MulTi-Raies) spectropolarimeter in He D3 line 
allowing direct estimates of the magnetic field strength in the prominence material. Small 
scale horizontal, field aligned features were observed in the high-resolution SOT images, 
showing evidence of quasi-period oscillations and upward (radial) propagation in the foot. 


Schmieder et al. (2013) found that the propagating small-scale features move upward with 


velocities of 10 km s _1 , period ~300 s, and a wavelength ~2000 km, and interpreted the 
features in terms of linear fast magnetosonic waves. The interpretation was supported by 
estimates of possible ranges of various wave speeds in the prominence foot (sound, Alfven, 
fast magnetosonic) based on the observed field strength, estimated density and line of sight, 
and by linear 2D MHD model of trapped linear fast magnetosonic waves. However, the 
amplitude of the oscillations is not small compared to the background intensity and exhibits 
increased relative magnitude with height. The waves appear to exhibit sharp fronts, in 
disagreement with the linear model predictions of sinusoidal time dependence. This suggests 
that nonlinear effects are important in the dynamics of these waves, affecting steepening 
that in turn can result in faster dissipation of the energy flux carried by the waves. 


In the present study we extend the work of Schmieder et al. (2013) by performing 


additional detailed analysis of the 2012 Oct 10 prominence focusing on the observed 
propagating wave-like features at lower heights in the prominence foot, and finding further 
evidence of their nonlinear nature. The results of our analysis are supported by nonlinear 
2.5D (two spatial dimensions and three components of velocity and magnetic field) MHD 
modeling of the propagating, trapped, nonlinear fast magnetosonic waves. We find that the 
interpretation of the observed propagating features in terms of nonlinear fast magnetosonic 
waves is in better agreement with observations then previous linear studies. We also study 
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the nonlinear waves produced by the magneto-acoustic gravity inodes in a prominence. 
The paper is organized as follows: in Section [2] we present the observational results of the 
2010 Oct 10 prominence foot oscillations, in Section [3] we present the MHD model, the 
initial state, and the boundary conditions used in the present study, Section [| is devoted to 
numerical results, and Section [5] to discussion and conclusions. 


2. Prominence foot observation on 2012 Oct 10 


In Figure [T| the context image of the prominence observed in EUV on 2012 Oct 10 at 
03UT by SDO/AIA is shown. Figure [l^, shows the 193A image that corresponds to <1MK 
plasma emission. The image shows the quiet Sun region on the solar disk, an active region 
in the top right end of the limb, and the overlaying EUV loops off-limb. Some prominence 
material is evident as dark absorption features in the central off-limb region, marked by the 
arrow. Figure [lja shows the same region as seen in 304A emission that corresponds to cooler 


material at ~50,000K emission (Labrosse & McGlinchey 2012). Here, the off-limb structure 
of the prominence becomes evident, and the location of the prominence foot is marked 
by the arrow. The animations in the online version of this figure show the various flows 
associated with the prominence material. However, the ~1” resolution of the SDO/AIA is 
not sufficient to show the propagating small-scale features discussed below. 


The Hinode/SOT image of the prominence foot in the Ca II H line on 2010 Oct 10 at 
14:04:47 UT is shown in Figure [2] The resolution of this image is about 0.1”, an order of 
magnitude better than SDO/AIA, showing clearly horizontal features with quasi-periodic 
separation. Part of the solar disk with spicules is evident as the saturated white region 
in the top-left corner. The online animation of this figure shows the dynamics of these 
features indicating that they are radially propagating wave-trains. In Figure |2jo a close-up 
of the prominence foot is shown. The black line indicates the location of the space-time 
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plots shown in Figure [3] below, and the saturated white region above the disk shows the 
emission from spicules. 

In Figure [3] the space-time plot of the emission along the line segment in Figure |2^i is 
shown. The heights are measured relative to the solar limb in the SOT. The propagating 
fluctuations are evident as bright and dark slanted features. The propagation velocities 
determined from the slopes are 8.5±1.2 km s^ 1 , 6.5±0.5 km s _1 , and 5.9±1.0 km s -1 for 
features 2, 5, and 6, respectively. The error bars of the velocities are estimated by the 
parallelogram method by fitting parallel straight lines that encompass the brightest features 
of the intensity perturbation. Since the temperature of the prominence material is on the 
order of 8,000K the sound speed is about 9 km s -1 . Since the magnetic field as estimated 
from THEMIS instrument observations is 5 — 14G in the prominence foot and the density is 
in the possible range of 10 9 — 10 12 cm -3 (with additional uncertainty due to the plane-of-the 
sky projection effects), the fast magnetosonic speed could be in a broad range of values (see 
Figure 11 in jSchmicder et al. 2013). Therefore, to narrow the estimates one needs to make 
assumptions about the parameters, such as the line-of-sight angle, and density that define 
the fast magnetosonic speed. For example, an estimate on the line-of-sight angle of the 
prominence can be obtained by comparison with realistic magnetic models of prominences 


e.g., Dudik et al. 2012). 


The temporal evolutions along the cuts denoted by the white lines in Figure [3] are 
shown in Figure [4] at the various heights in the prominence foot. The heights are at (a) 
11.1”, (b) 11.9”, (c) 12.7”, (d) 13.4”, (e) 14.2”, (f) 15.0”. The local peaks of the oscillations 
are numbered in each panel, and the non-sinusoidal and quasi-period structure of the 


intensity variations is evident. Schmieder et al. (2013) performed the wavelet analysis of 


similar oscillations and found periods of ~ 5 minutes. Here we are looking at variations at 
lower heights in the prominence foot as shown in this figure. 
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We have investigated the periods of the oscillations by performing wavelet (Torrence & 


Comp o|[T 998) and periodogram (Scargle||l982) analyses on the time sequence of the intensity 


fluctuations. In Figure [5] we show the wavelet analysis of the temporal evolution of the 
normalized intensity in the prominence foot in the Ca II H line of Figure [4] at heights shown 
in Figure |4j The solid curves show the boundaries of the cone-of-influence for the wavelet 
analysis located at ±y/P from the boundaries of the time interval. The line plots show 
the corresponding global wavelet for each position. It is evident that the oscillations are 
coherent over short periods of time compared to the overall time interval of the oscillations, 
and the global wavelet shows broad peaks. In order to narrow down the estimates of the 
periods we have also applied the periodogram analysis to the time sequences (Figure [6]) at 
the six locations. The peaks of the periodograms in Figure [6] that correspond to the broad 
peaks of the global wavelets shown in Figure [5] are at (a) 8.2±1.4, (b) 8.Oil.9, (c) 11.3T1.2, 
(d) 9.8T1.4, (e) 5.0T0.5, and (f) 4.9T0.5 min, respectively. The error bars of the periods 
are estimated from the half-width at half maximum of the periodograms. 


Since a compressional MHD wave is expected to steepen with height due to the 
gravitational stratification of the background density, and conservation of the wave energy 
flux, we plot the height dependence of 61 /I of the waves, where 61 is the perturbed intensity 
in Figure [7} It should be noted that while some gravitational stratification of the density is 
expected, it is diminished by the effects of the background magnetic held that supports the 
prominence material, and it may have nonuniform height dependence due to the complexity 
of the magnetic and density structure of the prominence. The error bars of 61/I are shown 
with the vertical segments. We kncl in the examined features that there is an increasing 
trend in 61 /I with height, which is qualitatively proportional to 6n/n in the observed Ca II 
emission. It is evident that feature 2 shows an increasing trend above 12”, and feature 3 
shows an increasing trend between 11” and 14”, while features 4 and 5 show increasing 
trends between 11” and 13.5”. Moreover, it is evident that 61/I is not <C 1 with nonlinear 





effects becoming more important as the ratio increases. Thus, the height evolution of the 
observed propagating relative fluctuations is in agreement with the expected evolution of a 
nonlinear fast magnetosonic wave propagating into gravitationally stratified atmosphere. 


3. MHD Model 

We solve the nonlinear resistive MHD equations in two spatial dimensions with the 
standard notation for the variables given by 


dp 

— + V ■ 
dt 

(P V ) 

= 0 , 



(1) 

dw , 

• V)v 


E u 

1 J x B 

( 2 ) 

« + < v 

— 

P 

Vp + 

r r r 2 p 

<9B „ 

~dt = 

X (v 

xb; 

) + 

5 _1 V 2 B, 

(3) 

( d 

U +v 

■ v ) 

^ b 3 

11 

= 0. 


(4) 


The normalization in equations (JTj) (|4]) is given by r — * r/R s , t —>■ £/ta, v —> v/Va , 

B —* B/B 0 , p —> p/po, and p —* p/po , where R s is the solar radius, ta = R s /Va is the 
Alfven time, Va = B 0 /y/4:7rp 0 is the Alfven speed, B 0 is the background magnetic field, p 0 
is the background density, and po is the pressure in the corona outside the prominence foot 
at r = 1 R s . In the present study viscosity is neglected, as well as the effects of heating 
and cooling. Other physical parameters are the Lundquist number S (in the present study 
we set S = 10 5 and the resistivity is negligible), the Froude number F r = V\R S /(GM S ), 
where G is the gravitational constant and M s is the solar mass, and the Euler number 
E u = PoUp^a) = C sh V A> where C s is the sound speed in the present study. In the 
present model we set B 0 = 5 G, n 0 = 10 9 cm -3 , T 0 = 10 6 K, which results in Va = 345 
km s _1 , ta = 33.8 min, and C s = 166 km s _1 with 7 = 5/3. With this normalization, the 
thermal to magnetic pressure ratio f3 = 8 ttp 0 /Bq = 2 E u = 0.2776. The equations are solved 
in two spatial dimensions keeping three components of the velocity and the magnetic field 
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(2.5D) using Cartesian geometry with the 4th order Runge-Kutta method in time, and 4th 
order spatial differencing on a 256 2 grid. A fourth order numerical viscosity is included 


e.g. Hamming 1973) for stability purposes. The numerical code used in the present study 


is adapted from the code initially developed to study waves in coronal holes and plumes 


(Ofrnan & Davila 1997, 1998 Ofman et al. 1999 Ofman & Deforest 2000). The nonlinear 


solutions are obtained with respect to an initial state and for the boundary conditions 
described below. 

In optically thick plasma the relation between the density and the intensity is complex 
and nonlinear, and in principle should be modeled by computationally costly radiative 


MHD codes (e.g., Gudiksen et al. 2011). However, radiative MHD models were not yet 
applied to this observation, and it is not the goal of the present study to determine the 
exact values of the density or temperature, but to study the wave properties which may 
provide qualitative agreement with the present, more simplified MHD model. In order to 
relate the results of the observations to the MHD model calculations we make the working 
assumption that SI/I ~ Snjn. The proportionality assumption of the ratios implies to first 
order in the relative perturbation that / ~ n“, and covers the possible cases of optically 
thin (a = 2 for collisional excitation) as well as qualitatively optically thick (0 < a < 2) 
plasma. Nevertheless, the exact values of a or n do not affect our results. The first order 
approximation is valid due to small values of 51 /I in most cases, evident in Figure [7J 


3.1. Initial State and Boundary Conditions 

3.1.1. Constant horizontal magnetic field 


Observations and models indicate that the magnetic field inside the prominence foot 


is dominated by the horizontal magnetic component (e.g. Dudfk et al. 2012). Therefore, 
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in order to model the nonlinear waves in the prominence foot we initialize the model 
with uniform horizontal magnetic field B = B 0 x and nonuniform temperature and density 
profiles along the field of the form 


Tt'A — T — (T — T ■ \ P -\.^ x - x o)/ w ] A 

± ± max y- 1 -max - 1 min )° 5 

(5) 

n{x) = P 0 /T(x), 

( 6 ) 


where T min is the minimal temperature inside the foot, T max is the temperature outside 
the foot, w = 0.1 is the half-width of the foot, and x$ = 0 is the location of the center of 
the foot. In this study we assume T min /T max = 0.01. The corresponding density profile 
is obtained from the pressure balance conditions in the x direction using the normalized 
equation of state with constant thermal pressure pq. Thus, in the prominence foot the 
sound and the Alfven speeds decrease by a factor of 10 compared to the outside region. 
The x-dependence of the density and temperature for the above initial state are shown 
in Figure [ 8 } This configuration is similar to the prominence model of JoarderAz Roberts 


(1992b), but with continuous temperature and density variation at the interface between 
internal and external regions. 


The fast mode waves are driven by periodic velocity perturbations at the lower 
(coronal) boundary by imposing time dependent fluctuations inside the foot given by 

V z (x, Z = 0,t) = y {cosut + l) e —[(*— *0 )/ 2 w}^ (7) 

where Vo is the amplitude of the velocity perturbation, u = 26.545 is the normalized 
driving frequency that corresponds to an 8 min period - close to the average period of the 
observed waves reported in Section [2] above. This form of the time dependence provides 
V z > 0 with Gaussian x-dependence centered at the model prominence foot. Symmetry 
boundary conditions are applied at x — 0 , with open boundary conditions at z = z rnax 
and at x = x max . The solutions are obtained in the 2D region (0,0.2) x (1,1.2) without 
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gravity. Thus, the gravitational steepening with height is not modeled. However, nonlinear 
dispersive effect and change in the background density of the waveguide are included. 
Hereafter, we refer to this model as ’model A’. 


3.1.2. Two dimensional magnetic field 






In the second part of our study we include gravity, using 
dimensional background magnetic field based on the prominence magnetic field model of 


Terradas et al. (2013) given by 


B x — BiCOskixe klZ — B 2 cosk 2 x e k2Z , ( 8 ) 

By = 0. (9) 

B z = — Bisinkix e~ kiz + B 2 smk 2 x e~ k2Z , (10) 


where k\ — k\ — |^, and L = 20. This form of the magnetic held provides B z (x = 0) = 0 
and the held is dominated by B x near x — 0 - the location of the model prominence foot 
(see, Figure [9]). The curvature of the held combined with gravity leads to trapping of the 
cool prominence material in the foot. Note, that for numerical convenience reasons the 
distances are normalized by 0.1/? s in this section, and the corresponding Alfven time is 3.38 
min. We also use here the polytropic energy equations with 7 = 1.05 to account implicitly 
for heating of the plasma column. 


The above initial state is supplemented by gravitationally stratihed density structure 
that in Cartesian geometry is approximated as 

n{z) = n 0 e“( r o+ 2 r o), ( 11 ) 


where n$ is the normalized density at z = 0 , ro = 10 is the location of the chromosphere- 
corona interface in units of 0 . 1 i? s , and a = GM s mn /( 2 /csToi? s ) is the normalized inverse 
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gravitational scale height, ks is Boltzmanns’ constant, rn h is the hydrogen mass, T 0 is the 
background temperature. The foot is initialized by introducing the temperature profile give 
by Equation [5] with T min /T max = 0.1 with the corresponding ^-dependence of the density 
multiplying the gravitationally stratified density, Equation 11 (see, Figure [13]) leading to 
over-dense prominence foot. Due to limitations of stability we use an order of magnitude 
lower temperature and density ratio than in model A. While the above background state is 
initially in equilibrium, the introduction of low temperature and high density region of the 
prominence foot results in gravitationally unstable initial configuration. This results in the 
initialization of the oscillations (instead of Equation [T] of model A). The solution is obtained 
in the 2D region (—2,2) x (0,4) with the following boundary conditions: line-tied magnetic 


field at z = 0 and open boundaries at the three other planes. While the model of Terradas 


et al. (2013) includes localized prominence material suspended in the model magnetic field, 


the present model is aimed at studying the prominence foot, and is significantly different in 
the distribution of mass and the initial state from the above model. Hereafter, we refer to 
the present model as ’model B’. 


Numerical Results 


4.1. Model A 


In Figure [TO] we show the density and the magnetic field structure of the model A 
prominence foot using the 2.5D MHD equations with injected waves amplitude V z q = 0.02 


at t — 1.68. Figure 10 (top) shows the density structure of the foot in the symmetric 
half-plane with magnetic field lines denoted by white lines. The density enhancements 
due to the propagating fast magnetosonic waves driven by the periodic velocity injections 
(Equation [7]) at the coronal base are shown. It is evident that the wavelength of the 
propagating density structures is ~ 20 Mm, and with the period of 8 min results in the 
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phase speed of about 42 km s^ 1 for the nonlinear waves, in good agreement with the 
theoretical value for linear transversely propagating fast magnetosonic wave inside the foot 
Vfi = (V^-t-C^) 0-5 = (34.5 2 + 16.6 2 ) 0,5 = 38.2 km s _1 . The values of Va% and C si are reduced 
by a factor of 10 compared to the values of Va and C s outside the foot due to density and 
temperature profiles shown in Figure [8] The corresponding magnetic field direction vectors 
and the magnetic field magnitude \B\ are shown in Figure [!o| (bottom). It is evident that 
the magnetic field enhancements are in phase with the density enhancement as expected for 
the fast magnetosonic wave. The wave pressure of the injected moderately nonlinear waves 
results in gradual modification of the background structure of the model prominence foot 
and consequent displacement of the magnetic flux as evident by the fieldline structure at 
the lower part of the foot. The propagating waves are mostly confined to the high-density, 
low-temperature foot, with some leakage evident in the magnetic structure. Animations of 
these figures are available in the online journal. 


The temporal evolution of the velocity components, perturbed density, n i, perturbed 
temperature, Tj, and the perturbed magnetic field components at height z = 1.1 in the 
center of the model prominence foot are shown in Figure [lT] for the driving velocity 


amplitude V z0 = 0.01, and at a height z — 1.1 for V z0 = 0.02 in Figure [12] The variations 
in the temperature perturbation are in-phase with the density perturbation as expected 
for the propagating wave. It is evident that the periodic velocity injection produces the 
fast magnetosonic wave with V z oscillations in phase with density, temperature, and B z 
oscillations. In addition nonlinear compression due to the waves produces V x that has a 
growing (secular) and oscillating part in anti-phase with the fast magnetosonic wave. The 
corresponding B x component is in quarter wavelength phase shift with the V x suggesting 
an Alfvenic nature of the nonlinearly driven secondary wave. The perturbed quantities are 
growing in time due to the nonlinear modification of the background structure of the foot 
due to the effects of the wave pressure. The y components of the velocity and and the 
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magnetic field remain zero in this 2.5D model. 

Increasing the amplitude of the injected waves by a factor of two (I4o = 0.02) makes 
the nonlinearity significantly more apparent in the variables. It is evident that the wave 
fronts steepen considerably with asymmetric fluctuations with shock-likc structures. The 
phase relations between the various velocity components are not affected and agree with 
the fast mode waves as in the less nonlinear case. The structure of the temporal evolution 
of the waves is in qualitative agreement with the temporal evolution of the intensity of 
the observed nonlinear fast magnetosonic waves shown in Figure [4] that show evidence 
of nonlinear steepening and non-sinusoidal shape of the waves at the various locations. 
The effects of the wave pressure on the background structure become more significant 
with increased wave amplitude. While 5(1) depends both on temperature and density 
oscillations, it is evident from the model that the variation of Tj is small, and the main 
temporal evolution of 5(1) is affected by the density oscillations. 


4.2. Model B 


In this model we show the results of the model B prominence foot where the initial 
magnetic field structure varies in two dimensions with dipped field as evident in Figures [9] 
and [13] The resolution in the x-direction in the 2.5D MHD model is doubled to 512 
since the full range of x is included in the model (i.e., the symmetry conditions are not 
applied at x — 0). In addition to the background equilibrium potential magnetic field and 
gravitationally stratified density, ‘heavy’ prominence foot material is introduced at t — 0 


as shown by the bright (high density) region in Figure 13. This results in gravitationally 


unstable initial state that produces gravity mode oscillations where the initially potential 
dipped magnetic field is dragged down by the heavy prominence foot material that is 
(nearly) frozen-in to the field in most of the structure due to low resistivity. This in turn, 
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leads to stretching and bending of the field and the formation of currents that produce a 
restoring Lorentz force, leading to oscillations of the magnetized foot structure (see online 


animation of the Figure 13). Since the lower part of the magnetic structure contains an 


x-point, due to the change of the magnetic topology, the prominence foot material can 
slip in this region due to the finite diffusion and gravitational acceleration affecting the 
oscillations, eventually destabilizing and disrupting the high density structure. 

The nonlinear gravitational mode oscillations are dramatically different from the 
nonlinear fast magnetosonic wave described above. In Figure [14] the oscillations at a 
height of z — 0.035i? s at the center of the prominence foot are shown. It is evident 
that the oscillations are dominated by the vertical component of the velocity 14 and the 
corresponding magnetic field component B z in the center of the foot, and that they are a 
quarter period out of phase. The density perturbation oscillations are in phase with V z 
with more complex dependence on the small Tf. The period of the oscillations is much 
longer (~6 hrs) than the fast magnetosonic waves, even accounting for the different size 
and parameters of model B with respect to model A. The shown duration of the oscillation 
is before the disruption of the heavy prominence material due to reconnection at the foot. 
Observational evidence suggest that qualitatively similar oscillations with gravity as the 


restoring force are observed by SDO/AIA (e.g., Li & Zhang 2012 Luna et al. 2014). 


5. Discussion and Conclusions 

Recent high spatial and temporal resolution observations of a prominence foot with 
Hinode/SOT Ca II emission show evidence of upward propagating disturbances. THEMIS 
observations provide the diagnostic of the prominence foot magnetic field and find fields in 
the range 5-14 G. The typical density of the prominence material in the range 10 10 -10 12 
cm -3 and temperature of ~8000K provide constraints on the possible sound, Alfven, and 
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fast magnetosonic speeds. These observations were interpreted as linear fast magnetosonic 


waves by Schmieder et al. (2013). 


We find that the disturbances are likely trapped nonlinear fast magnetosonic waves 
propagating perpendicular to the magnetic field of the prominence foot. This conclusion 
is supported by observations of strong nonlinear features, such as the non-sinusoidal form 
of the oscillations with sharp fronts, and an increase of SI/I ~ Sn/n with height over a 
significant range of heights in the prominence foot. The assumed proportionality of the 
ratios allows to first order a power law dependence between the observed intensity and the 
density of the pillar that may approximate both, optically thin and thick plasma. However, 
the exact values of the (positive) power or the density do not affect our results. 


We perform 2.5D MHD modeling of the nonlinear fast magnetosonic waves in a 
model prominence foot with horizontal magnetic field, high density, and low temperature 
(model A) using the typical parameters of the observations. The waves are driven by 
periodic velocity upflows at the coronal boundary producing associated density and 
magnetic field compressions. We find qualitative agreement of the modeled nonlinear fast 
magnetosonic waves features with observational signatures from Hinode/SOT that supports 
our interpretation in terms of nonlinear waves. The nonlinearity affects the propagation of 
the waves and the possible energy flux that can be carried and dissipated more rapidly than 
by linear waves. 


Using the 2.5D MHD model we also study nonlinear MHD waves excited in a model 
prominence foot 2D spatially variable magnetic field due to the effects of gravity on the 
heavy prominence material supported by the curved magnetic structure (model B) and find 
substantially different large scale and slower evolution from the nonlinear fast magnetosonic 
waves. The gravity-MHD waves in the model 2D magnetic field configuration are global 
(i.e., distant parts of the prominence foot oscillate in phase), exhibit much longer periods 





than the fast mode waves, and eventually result in the destabilization of the prominence 
material. These nonlinear waves are extensions of the well known linear magneto-acoustic 
gravity modes in prominences. Although, we have explored this possible wave mode 
numerically, we find that the the gravitational waves modeled here are not associated with 
the small-scale fluctuations seen by ffinode/SOT, but to global oscillations, reported in 
previous observations. 
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would like to acknowledge funding received through the NASA Earth and Space Science 
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Fig. 1.— The context image of the prominence from SDO/AIA in EUV on 2012 Oct 10 at 
03UT. (a) 193A , the arrow marks the location of absorption due to prominence material 
(b) 304A, the arrow shows the location of the prominence foot. Animations of this figure 
are available in the online journal. 





- 22 - 


-320 

-340 


S 
>■ 

-420 

880 900 920 940 960 980 

X (arcsecs) 


-360 

-380 

-400 


in 

o 

<u 

</> 

o 



900 910 920 930 940 

X (arcsecs) 



Fig. 2.— (a) The prominence in Ca II H line obtained with Hinode/SOT on 2012 Oct 10 
14:04:47 UT. (b) Closeup of the prominence foot (shown with black box in (a)) where the 
oscillations were observed. The black line shows the location of the space-time plot. The 
lowest part of the prominence is obscured behind the limb. Animations of this figure are 
available in the online journal. 
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Fig. 3.— Space-time plot of the prominence foot in Ca If H line obtained with Hinode/SOT 
on 2012 Oct 10 along the line shown in Figure |2Jd. The numbers refer to the locations of 
the peaks of the propagating waves and the white lines show the locations of the temporal 
evolution in Figure [4j The red dashed lines show several typical slopes of the propagating 
features, used to estimate the propagation velocity: 8.5±1.2 km s^ 1 , 6.5±0.5 km s -1 , and 
5.9±1.0 km s -1 for features 2, 5, and 6, respectively. The error bars of the velocities are 
estimated by the parallelogram method. 
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Fig. 4.— The temporal evolution of the normalized intensity in the prominence foot in the 
Ca II H line obtained with Hinode/SOT on 2012 Oct 10 along the lines shown in Figure [3] 
The temporal evolution at heights of (a) 11.1”, (b) 11.9”, (c) 12.7”, (d) 13.4”, (e) 14.2”, (f) 
15.0”. 
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(c) Cut at 12.7 arcsecs 
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(e) Cut at 14.2 arcsecs 
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Fig. 5.— The wavelet analysis of the temporal evolution of the normalized intensity in the 
prominence foot in the Ca 11 H line shown in Figure [4] at heights of (a) 11.1”, (b) 11.9”, (c) 
12.7”, (d) 13.4”, (e) 14.2”, (f) 15.0”. The curves show the boundaries of the cone-of-influence 
for the wavelet analysis (located at ±\/P from the boundaries of the time interval). The 
line plots show the corresponding global wavelet for each position. 
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Fig. 6.— The periodogram analysis of the temporal evolution of the normalized intensity in 
the prominence foot in the Ca II H line shown in Figure [4] at heights of (a) 11.1”, (b) 11.9”, 
(c) 12.7”, (d) 13.4”, (e) 14.2”, (f) 15.0”. The peaks are at periods 8.2±1.4, 8.0T1.9, 11.3T1.2, 
9.8±1.4, 5.0±0.5, and 4.9±0.5 min, respectively. The vertical dashed lines correspond to the 
peaks that are in agreement with the peaks of the global wavelets shown in Figure [5] 
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Fig. 7.— The height dependence of the relative intensity 51/I of the upward propagating 
disturbances for three typical features marked on Figure [3] in the prominence foot, (a) 
Feature 2, (b) feature 3, (c) feature 4, (d) feature 5. It is evident that features 2 show 
an increasing trend above 12”, feature 3 shows an increasing trend from 11” to 14”, while 
features 4 and 5 show and increasing trend between 11” and 13.5”. 
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Fig. 8.— The initial scaled density and temperature profiles across the magnetic field 
model prominence foot in the 2.5D MHD model. The density is in units of no = 10 9 
x is in units of 70Mm, and the temperature in units of T 0 = 10 6 K. 
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Fig. 9.— The initial structure of the 2D magnetic field of the prominence model. The arrows 
show the direction of the field, and the intensity shows the normalized magnitude. 
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Fig. 10.— Top: the density and magnetic field lines of the model prominence foot in the 2.5D 
MHD model at t — 61.8 min showing the structure of the propagating fast magnetosonic 
wave for the case with V z q = 0.02VSi = 6.9 km s -1 . The color bar shows the density range. 
Bottom: The corresponding magnetic field direction vectors and magnitude. The color bar 
shows the magnetic field magnitude range. Animations of these figure are available in the 
online journal. 
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Fig. 11.— The temporal evolution of the variables at height z = 0.1i? s at the center of 
the model foot. The initial transient solution at t < 34 min is not shown, (a) Velocity 
components, (b) perturbed density (ni), (c) perturbed temperature (T)), (d) perturbed 
magnetic field components. The velocity and magnetic field components are z (solid), y 
(dashes), x (dot-dashes). The driving wave amplitude is V z q = O.OIVa = 3.45 km s -1 . 
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Fig. 12.— Same as Figure 11 but with driving wave amplitude V z q 


0.0214 = 6.9 km s _1 . 
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Fig. 13.— The structure of the prominence foot density in the 2D magnetic held of the 
prominence model for T min /T max = 0.1 at t = 889 min. The white lines show representative 
magnetic held lines. Animation of this hgure is available in the online journal. 
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Fig. 14.— The temporal evolution of the variables (velocities, perturbed density, tempera¬ 
ture, and magnetic field components) at a height z = 0.035i? s at the center of the promi¬ 
nence foot for the case shown in Figure [l3j The velocity and magnetic field components are 
z (solid), y (dashes), x (dot-dashes). 








